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Abstract. During Run II of the Tevatron collider, D0 collaboration has made extensive searches 
for the neutral MSSM Higgs bosons (0), produced in pp collisions at -y/s = 1.96 TeV. Two such 
analyses, addressing inclusive <f> production with — > t + t~, and associated <j>b(b) production with 
<p — > are reported here. No excess of events above the background expectation has been observed 
in any of these analyses. The results are combined to set constraints on the MSSM parameter space. 
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: INTRODUCTION 
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q 1 The Standard Model (SM) of particle physics has been observed to be a consistent 

theory of fundamental particles and their interactions up to the energies they have 
been studied although one remaining essential ingredient the Higgs boson is yet to be 
discovered. Theoretical considerations teach us that SM can not be the ultimate theory of 
elementary particles and their interactions. Amongst many other compelling theories, the 
Minimal Supersymmetric extension of the SM (MSSM) is one of the viable theoretical 
frameworks which has the potential to overcome some of the short comings of the SM 
■ by incorporating minimal supersymmetric particle spectrum. 

Because of its supersymmetric structure, the MSSM[[l|][|2] requires at least two Higgs 
doublets to generate masses to both "up"-type and "down"-type quarks (and the re- 
spective charged leptons). Such a two-Higgs-doublet model predicts five physical Higgs 
bosons: two CP-even Higgs bosons, h and H, one CP-odd Higgs bosons, A and a pair 
of charged Higgs bosons, H ± . The MSSM Higgs sector at the lowest order can fully 
be described in terms of one Higgs boson mass (Ma is chosen in CP-conserving sce- 
nario and M H ± in CP- violating scenario) and tan/3 = V2/V1, where V2(vi) refers to the 
Higgs field that couples to the "up"("down")-type quarks. However, additional param- 
eters viz., Msusy, M2, ll and trig enter at the level of radiative corrections. It is to be 
noted that the coupling of the Higgs boson A to "down"-type quark such as bottom quark 
is proportional to tan/3. So with respect to the SM, the production cross section for A 
in association with bottom quark(s) gets enhanced by a factor of tan 2 /3. Furthermore, at 
large tan/3 (~ 50), there is mass degeneracy between the Higgs bosons A and h or H 
depending on the value of dia. Therefore, the total production cross section for MSSM 
neutral Higgs bosons i.e., h/H/A = is twice that of the Higgs boson A. Previously, 
LEP experiments have set the mass of light Higgs boson mh to be greater than 92.8 
GeV[3] at95%CL. 



The MSSM neutral Higgs bosons (0) mostly decay into bb pairs (90%) or into 
T + T~ pairs (~ 8%). So, at the Tevatron pp — ► 0(— > bb)b(b)X and — > 0(— > t + t~)X 
processes are considered to be the most promising channels to look for the signature of 
MSSM neutral Higgs bosons. Here we report on the searches for MSSM neutral Higgs 
bosons in the above channels performed by the D0 experiment [ 4] . The combined results 
are interpreted in different MSSM benchmark scenarios o]: "m™ ax " and "no-mixing". 

SEARCH FOR pp -» bb)b(b)X 

The analysis focuses on MSSM neutral Higgs boson production in association with one 
or two bottom quark(s) resulting in three or four bottom quarks in the final state while, 
CP-conservation in the Higgs sector is assumed. 260 pb 1 of D0 data collected between 
November 2002 and June 2004 have been utilized for this analysis. The data are first 
filtered using a dedicated on line trigger designed to maximize the signal acceptance. 
The Secondary Vertex (SV) tagging algorithm has been used to select the b-jets. The 
D0 SV b-tagging algorithm selects central b-jets (px > 35 GeV) with an efficiency of 
~ 55% while the mis-tag rate for similar light quark jets is about 1%. 

For Higgs masses between 90 and 150 GeV, the hb — > bbb signal events have been 
generated using the Pythia |6] event generator followed by the full D0 detector sim- 
ulation and reconstruction. The Pythia generated px and rapidity spectra of the Higgs 
bosons have been adjusted to those from NLO calculations [7]. The largest background 
contribution from QCD multijet processes have been determined directly from data 
while contributions from other background processes such as tt, Z°(^ M>)+jets have 
been determined from Monte Carlo. 

The selected events are required to have at least three b-tagged jets (within \f]\ < 
2.5) with px > 35, 20 and 15 GeV respectively for the leading px jets. The invariant 
mass distribution of the two leading px jets in data is then compared with that of the 
background expectation and no excess has been observed. Fig. Ufa) shows the limit 
on the production cross section at 95% CL, as a function of mA for tan/3 = 80 in "no 
mixing" scenarios. The results are also interpreted in the "maximal mixing" case and the 
limits in tan/3 — mA plane for both scenarios are shown in Fig.[TJb). 
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FIGURE 1. (a) The upper limits on the production cross section as a function of Higgs mass and (b) 
the limits on tanj3 — mA plane at 95% CL. 
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The analysis is based on 325 pb 1 of data collected by the D0 experiment during the 
period September 2002 to August 2004. The signal consists of a pair of tau leptons. 
One of the tau leptons is required to decay leptonically into electron or muon, leading to 
three final state signature: eXh, \iXh and efi, where T/, represents a hadronically decayed 
tau lepton. The three types of hadronically decayed x leptons, 

• T-type 1 : a single track with a calorimeter cluster without any electromagnetic sub 
clusters (1 -prong, 7T-like) 

• T-type 2: a single track associated with a calorimeter cluster and electromagnetic 
sub clusters (1 -prong, p-like) 

• T-type 3: two or at least three tracks with invariant mass below 1.1 or 1.7 GeV 
respectively (3-prong) 

have been utilized in this analysis. A set of neural networks 1 , one for each type have 
been used for further discrimination over the backgrounds. 

The signal as well as various background processes have been generated using Pythia. 
Events are then passed through the full chain of D0 detector simulation and reconstruc- 
tion software. Apart from the QCD multijet background where a jet mimics an elec- 
tron/muon from tau lepton, all other background processes are normalized to NLO and 
NNLO (Z° boson, W ± boson, Drell-Yan and di-boson) cross sections. 

The reconstructed visible mass defined as M v i s = \f (PtI + Pt2 + Pt) 2 > where P T i 2 is 
the four vector of the visible tau decay products and Pp = (ET,Ex,Ey,0), is used to dis- 
criminate between signal and background processes. No excess over the background ex- 
pectation has been observed. The Figs.|2£a)-(b) show the distributions of M v i s (denoted 
as "Inv. Mass(e//i,tauhad 5 Sr)" and "Inv. Mass(e,/i,Er)" respectively in the figure) for 
different final states, while Fig.|2c) displays the upper limits on production cross section 
as a function of Higgs mass. 
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FIGURE 2. The distribution of reconstructed visible mass M v ; s (see text) for (a) eT/,+/XT/, channel, and 
(b) ejj, channel, (c) The upper limits on the production cross section as a function of Higgs mass at 95% 
CL. 



Identical to the ones used for Z / y* — ► T + T cross section measurements! 8] 



COMBINED LIMITS 

The D0 results obtained from the above two analyses Il9ll lll0ll are combined together 
and are re-interpreted in "m™ ax " and "no-mixing" scenarios. Fig. |3] shows the current 
D0 limits on MSSM parameter space. The similar limits obtained by the CDF 111 ill are 
also displayed for the purpose of comparison. 
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FIGURE 3. The excluded region in the (Ma, tan/3) plane for the m™ ax and no-mixing scenario with 
U I 200 CieV. 



CONCLUSIONS 

D0 has performed searches for the MSSM neutral Higgs bosons and the current D0 
limit on MSSM parameter space is the most sensitive to date. Additional search channels 
like pp — > 0(— > % + T~)b(b)X, advanced analysis techniques and larger dataset will 
provide further scope to improve the sensitivity in future. 
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